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In the last three decades, researchers have tried to estabhsh universal patterns about the structure 
of food webs. Recently was proposed that the exponent 77 characterizing the efficiency of the energy 
transportation of the food web had a universal value (77 = 1.13). Here we establish a lower bound 
and an upper one for this exponent in a general spanning tree with the number of trophic species 
and the trophic levels fixed. When the number of species is large the lower and upper bounds are 
equal to 1, implying that the result rj = 1.13 is due to finite size effects. We also evaluate analytically 
and numerically the exponent rj for hierarchical and random networks. In all cases the exponent r; 
depends on the number of trophic species K and when K is large we have that 77 — » 1. Moreover, 
this result holds for any number M of trophic levels. This means that food webs are very efficient 
resource transportation systems. 

PACS numbers; 87.10.+e, 87.23.-n 



Understanding energy and material fluxes through 
ecosystems is central to many questions in ecology 

Ecological communities can be studied via resource 
transfer in food webs . These webs are diagrams show- 
ing the predation relationships among species in a com- 
munity. Usually, a group of species sharing the same 
set of predators and the same set of prey is aggregated 
in one trophic specie |E @- So, each trophic specie is 
represented by a site and denoted by an integer number 
i = 1,... ,K, where K is the total number of trophic 
species. A relation between a pair of sites is represented 
by a link directed from the prey to predator. There are 
several quantities introduced in the literature to char- 
acterize the food web structure, such as the fractions of 
the species in the trophic levels (basal, intermediates and 
top), the fractions of links among them, the connectance, 
the average distance between two sites, the clustering co- 
efficient and the degree distribution. It turns out that 
all these quantities are nonuniversal and dependent 
of the size of the food web. Perhaps, the only variable 
in common agreement with the literature is the maxi- 
mum value of trophic levels (M < 4). Garlaschelli et 
al. 1^ have considered food web as transportation net- 
works whose function is to deliver resources from 
the environment to every specie in the network. In this 
case, food webs appear to be very similar to other sys- 
tems with analogous function, such as river and vascular 
networks. In that work they have represented a real food 
web by spanning trees with minimal lengths. For each 
specie i the number Ai of species feeding directly or indi- 
rectly on i, plus itself, is computed. They also computed 
the cost of this transfer, namely C; = X^fe^fc' where k 
runs over the set of direct an indirect predators of i plus 
itself. The shape of Ci as a function of Ai follows a power 
law relation C {A) ~ A^, where the scaling exponent 77 
quantifies the degree of optimization of the transporta- 
tion network. They found the same allometric scaling 



relation for different food webs. By plotting Ci versus Ai 
for each one of the seven large food webs in the litera- 
ture, and by plotting Co versus ^0 for a set of different 
food webs. The exponent found, varies between 1.13 and 
1.16. Therefore, they concluded that the exponent 77 has 
a universal value (77 = 1.13) and it is, perhaps, the only 
universal quantity in food webs. Nevertheless, this mat- 
ter has been the subject of debates [iTlll^ . 

Here we establish an upper bound (r]max) and a lower 
one {rjmin) for the exponent 77 in a general spanning tree 
with M trophic levels and K trophic species, both fixed. 
In the limit iiT — > 00, we have that rjmax — Vmin — *■ 1- 
We also evaluate analytically and numerically the expo- 
nent 77 for hierarchical and random networks. Our main 
conclusions are that (a) the result rj = 1.13 for food webs 
is due to finite size effect (small K), (b) the exponent 
77 depends on the K and when K is large we have that 
77 = 1. Moreover, this results hold for any number M 
of trophic levels, implying that food webs are efficient 
resource transportation networks. 

It is worth mentioning that this problem is related 
to river and vascular networks 9]. Consider K sites 
uniformly distributed in a d-dimensional volume. The 
network is constructed by linking the sites, in such way 
that there is at least a path connecting each site to the 
source (a central site). Since each site is feed at steady 
rate Fi — F,the metabolic rate B clearly is given by 
B = J2i Pi = FK. Let lb represent the magnitude of flow 
on the h^^ link. Then, the total quantity of nutrients in 
the network, at a particular time, is given hy V — h- 
They deflne the most efficient class of network as that for 
which V is small as possible. Using this procedure they 
found that V ~ For river basins, d ~ 2 and 

V ~ B'^/^. In vascular systems V ~ B^^^ since d = 3. 
The variables Aq and Cq of the food webs are related, 
respectively, to the number of transfer sites N and to the 
total volume of nutrients V by the following equations: 
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N = Aq - 1 and V = Cq - Aq. Then we have that 
Co « A'^^^^'^^ if ^0 is large enough. The value of the 
exponent rj for a food web can be smaller than the one of 
rivers (ry = 3/2) or the one of vascular systems [t] — 4/3) 
because the spanning tree of a food web is not embedded 
in an Euclidean space. 

Let us consider a hierarchical network with M trophic 
levels. The network is constructed in the following way. 
We begin with a site representing the environment, the 
site 0. Then we connect ni sites to it, since these sites 
are feeding directly of the environment they constitute 
the first trophic level. Obviously, the number of species 
in this level is A^i = ni. The second level is constructed 
by connecting n2 sites to each site of the first level. Now, 
in this level, we have N2 = nin2 species. This procedure 
is repeated until the level M. 

Since Ai is the number of species feeding directly or 
indirectly on site i, plus itself, we have that 

Am = 1 
Am-1 = nMAM + 1 = TT-M + 1 

Am-2 = Um-iAm-I + 1 = nM"M-l + «M-1 + 1 
M 

Ao = l + ^7V^ = if+l . 

The cost of resource transfer, defined by = A^, 
where k runs over the set of direct and indirect predators 
of i plus itself is given by 

Cm = 1 

Cm-1 — ^mCm + = 2nM + 1 

Cm-2 = nM-iCM-i + Am-2 = SnMriM-i + ^um-i + 1 

Ma M 
a=l i=l (5=1 

The exponent 77, as was proposed in the literature @, 
can be found by (a) plotting Ci as a function of Ai for a 
network with number of trophic level M and total specie 
number K fixed; Usually, the point (1,1) is neglected 
due to finite size effects. It can be also found by (b) 
plotting Co as a function of ^o for several networks with 
different trophic species number K. This last procedure 
determines the large scale exponent |T2|. Note that in 
hierarchical spanning tree networks, Ci and Ai for species 
in the same trophic level are equal, implying that we 
have only M + 1 points in a C^ x At plot. Let us firstly 
use procedure (a) for networks with constant ramification 
ratio Ui = n and constant number of trophic levels M = 
4. We find r] = 1.39 for n = 2 and = 30 and = 1.27 



for n = 3 and ii" = 120, as it is shown in figure ^ Clearly, 
the exponent 77 depends on value of K, and decreases as 
long as K grows. In the limit that n — )■ 00, the total 
number of species K also is unlimited and we obtain that 
the exponent 77 approaches the value 1. 
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FIG. 1: Log-log plots of d versus Ai for networks with con- 
stant ramification ratio n — 2 and n — 3. Note that the 
exponent decreases when n grows. 

Let us return to the more general case of hierarchical 
models. The large scale exponent 77 can be evaluated by, 

InCo _ ln[l + Ef^i(^ + l)nti»d ^^ 

inAo in(i + E'ii nr=i «0 

If at least a ramification ratio is large, Uf} — > 00, we 
have that InAo ~ In 71^3 and InCo ~ Innp. Therefore we 
find r; — > 1 when the number of species is large. We can 
also use Eq.^to evaluate the exponent 77 for hierarchical 
networks with constant ramification ration. We find for 
this networks 77 = 1.41 {n = 2 and K = 30) and 77 — 1.31 
{n — 3 and K = 120). These values can be compared 
with the ones obtained previously with procedure (a) (see 
Fig.lIJ. 

In the Eq. ^ the exponent 77 depends on value of K, 
decreasing as long as K grows. For example, consider the 
hypothetical food web with total specie trophic number 
K — 146 and the specie trophic numbers in each level 
given by iVi = 38, N2 63, A^3 = 43 and = 2. We 
find the exponent r/ = 1.22. But, if we double the number 
of trophic species in each trophic level Ni = 2Ni the 
exponent is now 77 « 1.19. In that equation the exponent 
77 also depends in the relative distribution of the species 
in each level, for a given total specie number K. For the 
hypothetical food web described above with 146 trophic 
species we change the distributions of species in each level 
to iVi = 114, N2 = 20, N3 = 10 and A^4 = 2. We find 
the exponent 77 — 1.16. The exponent has changed from 
77 = 1.22 to 77 = 1.16 

Now, let us consider a random network with M trophic 
levels and K trophic species. The network is constructed 
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in the following way. First, we determine randomly the 
population in each level (a = 1, 2, . . . , M), obey- 
ing the restrictions M fixed and K fixed. Then, the iVi 
sites are connected to the environment, constituting the 
first trophic level. The second level is constructed by 
randomly connecting the N2 sites to the A^i sites of the 
first level. This procedure is repeated until the level M 
is constructed. In this case, we can evaluate the mean 
value of Ai and Ci in each level, namely 
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• Data for K= 123, M=4 
- Best fit, ri=1.18 
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Here a specify the trophic level (a = 1, . . . , Af). These 
quantities are averaged on several random configurations. 
Note that in the last level we have that Am ~ Ai = d = 
Cm = 1 and that we always neglect the point (1,1) in 
all fits 

In Fig. |2Ia) it is shown the Ca x Aa graph for a 
random network with K = 123, the same number of 
trophic specie that the Ythan Estuary with parasites, and 
M = A. A best fit furnishes = 1.18. A similar fit for 
K = 93, the same number of trophic specie that the Little 
Rock Lake food web, and M = 4 give us rj = 1.21. Note 
that the exponent decrease when K increases. Clearly, 
our exponent is larger than that found by ^ for the 
same trophic species number K. But, when K grows 
our exponent become smaller that them. Obviously, if 
rj = 1.13 represent a universal value for food webs of all 
the sizes, then random spanning trees networks with the 
same number of trophic levels M are more efficient than 
food webs. In fig. E^b) it is shown the Ca x Aa graph 
for a random network with K — 10000 and M = 4. Note 
that when K is large enough the exponent 77 « 1 

The exponent 77 can also be computed by the proce- 
dure (b) . For each value of K we perform an average for 
several configurations and find the mean value of Cq. In 
Fig. Ela), it is shown the Cq x Aq plot for random net- 
works with M — A and K varying from 50 up to 1000. 
Now we have that i] = 1.00. It is worth mentioning, 
that Cq X ^0 always furnishes 77 = 1 independently of the 
range of K. We have also simulated random networks 
with M = 10 trophic levels. In Fig. |2Ib) it is shown the 
Co X ^0 plot. The results are similar. 

Now let us present the central point of this paper, a 
general argument to demonstrate that the large scale ex- 
ponent is 77 = 1 for large K . Let us consider a spanning 
tree with M and K, both fixed. To obey the constraint 
of M fixed, we put one site in each level. Now we must 
put each one of the reminder K — AI sites. Since Cq is 
cumulative, a site put as near as possible of the envi- 
ronment has the minimal contribution to the global cost. 
On the other hand, a site put as far as possible of the 
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• Data for K= 1 0000, M=4 
— Best fit, ri=1.10 
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FIG. 2: Log-log Plots for random networks, (a) Ca x Aa for 
a network with K — 123 and Af = 4 {h)Ca x Aa graph for 
a random network with K — 10000 and Af = 4. Note that 
when K is large enough the exponent 77 « 1 

environment has a maximal contribution to Cq. To con- 
struct the network with maximum value of Cq, CQ^rnax, 
we must link all K — AI sites to the site of the last level. 
In this network we have Ni ^ N2 = ■ ■ ■ = Nm-i = 1 
and Nm = K — M + 1. CQ^min is obtained by linking 
the K — M sites directly to the site representing the en- 
vironment. In this case, we have that Ni = K — M + 1 
and N2 — = . . . = Nm = 1- Note that these con- 
structions are the closest networks to the star-like and 
chain-like ones, respectively, that obey the constraints of 
M and K fixed. Using the Eq. QJwe have. 



Co,™„ = l + 2K+—{M -I) 

M 

CQ^max = l + K(M+l) + —{l-M) 

Then, the lower and the upper bounds for the exponent 
77 are 

„ lnCo,„a. ln[l + K{M + 1) + f (1 - M)] 

Tlraa. - ^^^^ ^ " \n{K + I) 
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• Data for M=4 
— Best fit, ri=1.00 



When K ^ oo, we have that r]„ 
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• Data for M= 10 
- Bestfit, ri=1.00 
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FIG. 3: Log-log Plots for random networks, (a) Co x Aq for 
M = 4 and K varying from 50 up to 1000 and (b) Co x Aq 
for M = 10 and K varying from 200 up to 5000. Note that 
the large scale exponent has the same value rj — 1.00 and is 
independent of M. 



InCo^^.n ^ Hl + 2K+f{M-l)] 
ln{K + 1) \n{K + 1) 



Consider again the simulation of random networks. We 
verified that the constructions with minimum and maxi- 
mum Co are the ones just described. Moreover, the resuh 
above explains why we find 77 — > 1 when K is large in the 
simulations of random networks. 

In summary, we studied the transportation properties 
of several networks that represent spanning trees of food 
webs. First, we analyzing an idealized hierarchical model 
that can be analytically solved. Then we show that the 
exponent rj depends on value of K and, in the limit that 
K is large enough, the exponent 77 approaches the value 
1. After, we construct random networks that more real- 
istically represents a spanning tree formed by food webs. 
We evaluate numerically the exponent rj by several proce- 
dures. Again, in all cases the exponent depends on size 
of web and if K is large rj 1. One important point 
is that all the results are independent of the number of 
trophic levels M. Moreover, we establish a maximum 
and a minimum values for the exponent r/ in a general 
spanning tree with K and M fixed. When the number of 
species is large these values became equal to 1. There- 
fore, we must find rj = 1 for a large food web and we 
can conclude that food webs are very efiicient resource 
transportation systems. 
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